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Abstract: We give a local integral formula, valid on general curved space-times, for 
the characteristic Cauchy problem for the Dirac equation with arbitrary spin using the 
method developed by Friedlander in [S]. The results obtained by Penrose in the flat 
case in [50] are recovered directly. It is expected that this method can be used to obtain 
sharp estimates for the characteristic Cauchy problem for the Dirac equation. 
Resume: Nous donnons une formule integrale pour le probleme de Cauchy caracteris- 
tique local pour I'equation de Dirac pour le spin arbitraire en utilisant la methode 
developpee par Friedlander dans [H]. Nous retrouvons alors directement le resultat de 
Penrose dans le cas plat ([20j). 



Penrose obtained in 1963 an integral formula for the characteristic Cauchy problem for the 

Dirac equations for arbitrary spin in the flat case. His derivation of the integral formula is based on 
the construction of a Ncwman-Penrose tetrad (null tetrad) adapted to the null structure of the null 
initial data hypersurface, and especially to the description of the behavior of the null generators 
(bicharacteristics) of the cone. The use of the 2-spinor formalism allows him to write the solution 
of the problem in function of a "null datum", contraction of the data on the cone with its spinor 
generators. The formula is verified a posteriori through a splitting of the Dirac operator over the 
spin basis in the compacted spin coefhcient formalism. Penrose expected that this formula could 
be extended to the analytic case. As far as the author knows, the general case remains open. 

Friedlander gave in the mid 70's ([5]) a method to obtain a parametrix for the wave equation 
derived from the Leray constructions (see for instance [8 ) and wrote an integral representation of the 
solution of the characteristic Cauchy problem. His construction is based on a natural decomposition 
of the fundamental solution on the cone. Another approach exists to the characteristic Cauchy 
problem based on Fourier Integral operators. It must furthermore be noticed that there is no 
general result about the characteristic Cauchy problem for hyperbolic operators. Hormander gave 
in [13] a general result of existence and uniqueness, together with energy estimates, for the wave 
equation on a spatially compact Lorentzian manifold. 

The purpose of this paper is to combine the method developed by Friedlander with the descrip- 
tion of the null cone by Penrose to obtain an integral formula for the characteristic Cauchy problem 
with initial data on the cone for arbitrary spin in general curved spacetimes. The choice of this 
method implies that we face the same restrictions as in the book by Friedlander. There exists an 
essential obstacle to the extension of the domain of validity of the representation formula: the exis- 
tence of caustics which limit the domain where the formula can be written. We then have to restrict 
ourselves to a geodesically convex domain O of a smooth Lorentzian manifold (M, g), that is to say 
a domain where there exists a unique geodesic between any pair of distinct points. This restriction 
is inherent to the method and the fact we work with arbitrary curved geometry. The advantage 
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is however that we obtain an exphcit integral formula without resorting to any microlocalization. 
This in principle should allow an extension to metrics of low regularity in the spirit of [14j . 

More explicitly, let us consider {M,g) a smooth Lorentzian manifold and po a point in f2; the 
problem: 

= 

where m is a section of a given fiber bundle on and p is the Dirac operator on this bundle, with 
the initial conditions on the future null cone C+(po): 

u = 9 on e+(po) n Vt 

is known as a first order Goursat problem with initial data on the characteristic hypersurface 

e+(po)ni]. 

It is known that several conditions must be satisfied to ensure that this problem admits a 
solution. The first one comes from a geometric obstruction to the existence of a solution when 
symmetry conditions on the field u are imposed; this implies that the manifold M must satisfy 
some geometric assumptions, known as the consistency conditions, depending on the spin we are 
working with. The second one comes from the fact that the initial data are given on a characteristic 
hypersurface: 6 must then satisfy the restriction of the Dirac equation to the cone from pq: 

Ple+(po)^-0. 

These equations are called the compatibility equations for the initial data. 

As already mentioned, there exists, as far as the author knows, no general result about the 
characteristic Cauchy problem. Nonetheless, it is worth mentioning some results of existence and 
uniqueness with some generality. In the analytical case, this problem is similar to the Cauchy- 
Kowaleski problem (see for instance [8 ). The problem is well posed in that case. This can be 
extended, with energy estimates, to minimal regularity (|12j). The well-posedness of the character- 
istic Cauchy problem is nonetheless not the point of this paper: assuming existence and uniqueness 
of the solution in the neighborhood of the point the goal consists in deriving a representation 
formula for this solution. 

The paper is organised as follows. The first part presents an adaptation of the Friedlander 
method to the bundle of Dirac spinors. After a geometric and intrinsic presentation of the theory 
of spinors, the analytic tools to write a fundamental solution of the Dirac equation are developed. 
The second part is devoted to the derivation of the formula for Dirac spinors. Following Penrose's 
construction, a null tetrad adapted to the structure of the null cone is constructed and used to 
describe the geometric tools. The integral formula can then be derived from the parametrix and the 
result obtained by Penrose is recovered for Weyl (or two-) spinors. Finally, the third part deals with 
the arbitrary spin ^. The presentation made in the first part is adapted to the bundle of spinors 
with spin ^ so that the construction can be applied directly. A representation formula is then given 
for arbitrary spin and simplified in the case of the Maxwell equations. Penrose's formula for the 
characteristic Cauchy problem for arbitrary spin in the fiat case is recovered in a flat spacetime. 

The author would like to thank his supervisor, Jean-Philippe Nicolas, for his kind and patient 
support, his advices and his culture. 
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Notations and conventions. 

We describe here for future reference the notations and conventions which will be used all along 
the paper. Note that smooth means C°° in this paper. 

a. Geometric notations: 

(i) General framework: 

• {M,g): smooth Lorcntzian oriented and time oriented manifold with a metric g 
having signature signature (+,—,—,—); 

• fi: geodesically convex domain of M; 

• fj,: volume form associated with the metric g on M; 

• Pq: Si given point in Q; 

• V: Levi-civita connection for g on the tangent bundle of M, TM. 

(ii) Null structure on Q: let p be a given point in Q: 

• C{p): null cone from p, that is to say the set of points of il which lie on a null 
geodesic passing through p; 

• C^(p) (resp. C~(p)): future (resp. past) null cone from p, that is to say the set of 
points of fl that lie on a future (resp. past) oriented null geodesic from p; 

• 3(p): chronological set from p, that is to say the points of fl which lie on a timelike 
or null geodesic passing through p; 

• 3^(p) (resp. d^{p)): future (resp. past) chronological set from p, that is to say the 
set of points of Q that lie on a future (resp. past) oriented timelike or null geodesic 
from p; 

• 3ip) = 3{p)\e{p): causal set from p and 3^ip) = a=^(p)\e±(p) are the future and 
past causal sets from p. 

(iii) Spin structure: O is endowed with a spin structure; the spinors will be denoted using the 
Penrose conventions as well as the usual algebraic notations according to convenience: 

• ^Dirac- fibre bundle of Dirac (or 4-) spinors; 

• Sa and : bundles of Weyl (or 2-) spinors (resp. dual and anti-spinors) ; 

• "•": Clifford multiplication; 

• (•,•)■ symplectic product on ^Dirac obtained by lifting the metric g; 

• e^^ and eA'B'- restrictions of (•, •) to E>a and ; 

• C'^iSoirac) = "^{^Dirac)'- smooth scctious with compact support in Q endowed with 
the usual Frechet topology; 

• D' (EiDirac)- its topological dual; 

• C°°(SDiroc) = ^-i^Dirac)- smooth sections oiE>Dirac on fl,; 

• £'(§Dirac): its topological dual; 

• the connection V on TH. is lifted on §Dj>ac and is still denoted V; 
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• the Dirac operator is defined, for a given section (ei)ig{o,...,3} of the fibre bundle of 
orthonormal frames, on C°°(S£)irac) by: 

V$ e C°°{§D^rac), P$ = ^ ■ Ve. $ 

ie{0,...,3} 

1 Geometric and analytic preliminaries 

The geometric and analytic tools are presented in this section. As already mentioned in the intro- 
duction, due to geometric obstructions such as conjugates points or convergence of geodesies, the 
whole paper restricts itself to a geodesically convex domain Q: 

Definition 1.1 A domain Q is said to be geodesically convex if and only if it is an open set where, 
for every pair of points {p, q) in fl, there exists a unique geodesic between p and q. 

1.1 Dirac spinors and Dirac equation 

This section presents a construction of the spinor bundle so that it will be possible to apply the 
method of Friedlander in the most direct way. This presentation also intends to be a small dictionary 
between an abstract presentation of the theory of spinors and the Penrose conventions to represent 
spinors in terms of indices. Finally it must be noticed that, though the presentation is made on fi, 
it can be generalized to a globally hyperbolic manifold (see remark [T31 below) . 

1.1.1 Abstract construction 

We begin by defining a spin bundle: 

Definition 1.2 A manifold M is said to be spin if its tangent bundle admits a spin structure, that is 
to say there exists a Spin(l,3) principal bundle Ps, together with a twofold covering ^ : -Pg ^ PsoM , 
where PsoM is the SO {I, 3) -principle bundle of orthonormal frames on M, such that 

yp e Ps,yg e 5pz7i(l,3),e(ra) - 

where "i-s the universal covering from Spin(l,3) ~ SL2{C) on 50(1,3). 

Remark 1.3 a. The existence of a spin structure on a manifold is usually ensured by the as- 
sumption that its second Stiefel Whitney class vanishes. 

b. In the case of a four dimensional Lorentzian manifold {M,g), Geroch showed in \9I that 
a necessary and sufficient condition for M to carry a spin structure is that its bundle of 
orthonormal frames admits a global section (this is referred to as parallelizability) . 

c. A common assumption in general relativity which ensures that a 4^- dimensional Lorentzian 
manifold is spin is the global hyperbolicity assumption: there exists in M a global Cauchy hy- 
persurface, i.e. a spacelike hypersurface such that any inextendible timelike geodesic intersects 
this surfaces exactly once('9, 101). 

The spinor bundle on is defined through the action of an algebra over a vector space. This 
construction requires the following tool, which consists in group action over a fibre bundle, replacing 
its previous fibre by a given vector space: 

Definition 1.4 Let {E,VL,it) be a G-principal bundle. Let F be a vector space and 
p : G Homeo{F) a continuous map. 
Consider the action 

(j): G ^ Aut{E) X Homeo{F) 

9 ^ {{x,y) e E X F 1-^ {xg-'^,p{g)y)) 
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The quotient space 

E X F/p or EXpF 
with projection tt obtained by factorization of the diagram 

Ex F ^ M 





EXpF 

where p : E x F F is the projection on the first variable, is a G-principal bundle with fibre F . 

The algebra, known as the ChfFord algebra associated with a given quadratic form, which is 
used to construct the spinor bundle is then defined: 

Definition 1.5 Let E be a vector space (real or complex) with a quadratic form q. The Clifford 
bra {Cl{E), +, •) is the quotient space: 

/ +00 n \ 



0(g) 



(1.1) 



where I{E) is the ideal generated by the set {w (g) w — q{v)\v S E}. 
This algebra is known to have the following structure (|15j): 

Proposition 1.6 There exist two sub-algebrae denoted Cl^{E,q) and Cl^{E,q) such that: 

Cl{E, q) = Cl°{E, q) ® Cl\E, q) 

wich satisfy: 

Cl°{E, q) ■ Cf{E, q) = Cf{E, q), Cf{E, q) ■ Cl\E, q) = Cf{E, q) 
Cl°iE, q) ■ Cl\E, q) = Cl^ {E, q), Cl^ {E, q) ■ Cf{E, q) = Cf{E, q) 

Definition 1.7 The group Spin{E,q) is the subset of CP{E,q) defined by 

{seCf{E,q)\q{s)^\} 

where q is the extension of the quadratic form q to Cl{E,q). 

The formalism previously defined can of course be applied to the case of the Minkowski spacetime 

Definition 1.8 The bundle defined by: 

Smrac = (ft X M2{C))/{Spin{l,3)) 
is called the bundle of four dimensional spinors or bundle of Dirac spinors. 

Remark 1.9 a. The representation of Spin(l, 3) = SL2{C) acting on M2{C) has two irreducible 
components, which correspond to with its two inequivalent complex structures; by conven- 
tion, we write: 

M2{C) ^ (C^)* ® C2. 
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b. The previous remark gives a decompositon of the fibre bundle Soirac into two bundles (known 
as bundles of Weyl spinors), corresponding to the splitting of M2{C) into C and (C^)*/ this 
decomposition is written in terms of indices: 

^Dirac — © ^ 

A section u of the Soirac bundle will then be split into two smooth sections of the Weyl bundles: 

c. To get back to the tangent bundle, a convention must be chosen to represent the Clifford algebra 
Cl{M.^,ri). Its usual representation is M2(IHI), which is split in M2{C) ® A'/2(C). The vectors 
are identified with the hermitian 2-forms or with <C-antilinear homomorphisms from §a to 

. As such, a vector will be written by convention w^^ . 

d. For the chosen representation of the Clifford algebra which was made previously, the tangent 

bundle is identified to the set of hermitian two forms over S a- As such, it endows the tangent 
bundle with a structure of conformal Lorentzian manifold, i.e. a fibre bundle of cones over M 
and a time orientation: the fibre bundle of cones over SI is made of degenerate hermitian two- 
forms, the spacelike vectors fields are the hermitian matrices of signature (2,0) or (0,2) and 
the timelike vectors fields are the hermitian matrices of signature (1,1). The time orientation 
is obtained by a choice of orientation on the fibre bundle of cones over Q,. 

Proposition 1.10 The bundle A^S"* of skew- symmetric two forms is trivial. 

Proof This is a direct consequence of the fact that S"^ is a Spin{\,2>) — SL2{C) bundle. Let 
{U X C^, (p) and (V x C^,'^) be two local trivializations of the bundle S"^ with empty intersection 
where (j) and tp satisfy 

po (j) = n and poip = tt 

where tt : S"^' ^ M is the projection associated to the bundle S"^' and p is the projection on the 

first variable. These two trivializations give rise to two trivializations of A^S"^ that are still denoted 
by 4> and tp. Let us consider the transition map (potp^^ : A^C^ xV ^ A^C^ x U. It can be written: 

<P°ip{x,y) = {x,u{x)y) 

where v : U r\V ^ SL2{C) is a smooth map. 

Let X be fixed in U (IV. Since v{x) belongs to 5^2 (C) and y is a skew-symmetric 2-form, y is 
invariant under the action of an element of SL2{C), i.e: 

V(u, v) e C^, y{i>{x)u, i'{x)v) = y{u, v). 

The fibre bundle A^S"^' is thus trivial.® 

RemEirk 1.11 a. The canonical isomorphism, which will be denoted by k, between S"^ and §a 
induces an other isomorphism between A^S"^ and A^S^' 

A^S^' — > Is?^A 

e I — > /t,e : (u, ti) € X §^ e(K(w), 

It allows the construction of a symplectic form on §Diroc- l^t e in A^S"^ , we obtain a symplectic 
form on ^oirac by taking: 

e®K*e. 

A two-form e on is denoted s"^^ and acts on Weyl spinors by: 

y{UA, Vb) e Sa, v) = e^^UAVB- 

The corresponding two-form on §^ is denoted ea'B' ■ 
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b. Let £ he a fixed skew- symmetric two-form on §a- It is possible to construct a metric g on TQ 
by, for x^^ and y^^ two vectors: 

QabU^V^ = eABSA'B'X^^'y'^^' 

c. We denote by eab a two form which gives rise to the metric g on M. The non- degeneracy of 
£ induces an indentification between Sa and its dual S"^ given by: 

ka&Sa> — > = e^^KB e §^ 

whose inverse mapping is 

e §^ I > Kb = K^SAB- 

The equivalent transformation can be made for the complex conjugate spinors in if we 
consider the image two-form £a' b' ■ 

d. The symplectic product on Dirac spinors can thus be written, by lowering and raising indices: 

{u,v) = E^^'Mb+Sa'B^^'C''' 
= -^^A+U'C 

where u = tpA + and v = (pA + are two Dirac spinors. 

e. The dual ^^owac °/ ^Dirac is split in: 

The symplectic form (•,•) realizes an identification between ^oirac and S^jj^oc whereas its 

restrictions to, respectively, §>a and S'^ . denoted, e^^ and sa'B', realize an identification 
between and Sa and between Sa' and respectively. 

Proposition 1.12 Let e be a section o/A^S^. Let g be the metric associated with e. Then, the 

metric g is conformal to the metric g. 

Proof : Let p in M. Let X in TpM and u,v in We assume tliat tlie vector X is a liglit-lilce 
vector for the metric g, i.e.: 

g{X,X)=gai,X''X^ = Q. 

A necessary and sufficient condition for gab to be conformal to g is that g and g have the same null 
cone structure, i.e. it is sufficient to show: 

~g{X,X)=gabX'^X'' = Q. 

Since X°- is light-like vector, it can be written: 

if it is future directed and 

X'' = -u^u^' 

if it is past directed. The calculation is performed for a future directed null vectors, but it is the 
same for a past directed one. Because of the skew-symmetry of sab, we have: 

eabu'^u^ = 

and then 

gabX^X^ = eABU^V^EA'B'U^'u''' = 0. 

X is thus still a null vector for g and g and g are conformal metrics. 
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Remark 1.13 The map: 

A^S^ > {<j)g\(t>€C°°{n,Rl)} 

£AB I > gab = £AB£A'B' 

is a two sheeted covering of the conformal class of g. In particular, it is surjective. We denote by 
£ab a preimage of gab- 

Proposition 1.14 The bundle Soirac is a Dirac bundle, i.e. a fibre bundle of left modules over 
Cl{n,g) endowed with a symplectic form e and a connection such that: 

a. is the pull-back of the Levi-Civita connection on M : if tt : EiDirac — > then V can be 
written: 

7r*V = 

b. the connection is compatible with the action of the Clifford algebra: let X be a smooth section 
of Cl{Til, g) and u a smooth section ofE>Dirac, then: 

V^(X ■u)^yX-u + X- V^u. 

Though different since they are acting on different objects, the connexion V on f2 and on 
^Dirac o,i"e both denoted by V. 

c. the action of the Clifford multiplication is an isometry for the symplectic product: let X be a 
smooth section of Cl{Til, g) and u,v two smooth sections ofSDirac, then: 

e{X ■u,X-v) = q{X)e{u,v). 

In order to define the space on square integrable spinors onil, it is necessary to define the norm 
of a spinor. This unfortunately cannot be done without choosing a time function t (see |18j ) or, at 
least, a timelike vector field. 

Definition 1.15 A smooth function t on is called a time function if, and only if its gradient is 
a non-vanishing future- oriented timelike vector field on U. 

Definition 1.16 Let t be a time function on fl. Then the map defined by: 

(*,$) I y £(Vt-^',$) 

is a positive definitive hermitian product over the set of smooth sections of Soirac with compact 
support in U. The norm associated to this scalar product is denoted by \ \*\\u- 

Remark 1.17 • This norm will be used in the following on compact subsets of an open set of 
D, to define the Frechet topology over .smooth sections of the fiber bundle of Dirac spinors. 

• A time function t is fixed on fl. This time function will be used to compute all the norms. 

This scalar product is used to define various norms over the spinor fields on il: let <f> in Ti{SDirac)- 
We define using the positive definite hermitian product: 

• the L°°-norm over a compact K of fl: 



|$||oo,if = , /sup (£(Vt •$,$)); 
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if : O — > is a given chart over O, the norm over for any integer A'': 

||*||oo,JV,K = 



sup(£(Vt • V"<E>,V"$)), 



\ \a\<N 



where V" = Va^^ ^<5c<2 • • • ^^a, ' ~ ('^i' • • • ' '^i) being a multi-index of length 

1^1 = Si=i...n'^»' ^ chart of reference \1/ is fixed in the following in the computation of the 

norms; 

the i^-norm over f2: 



|^||2 = W^e(Vt-$,$)M. 



1.1.2 Newman- Penrose tetrad. 



One way to describe the Lorenztian structure is to use a global section of the fiber bundle of 

orthonormal frames over Q. and translate the result in terms of spinors. We construct a global 
basis, named tetrad of Newman- Penrose which gives rise to a spinor basis of S"^. 

Definition 1.18 A basis ofTVl®<C {l,n,m,rri) is called a normalized Newman Penrose basis if I 
and n are real vectors fields and it satisfies the following relations: 

g{l,l)=0 , 9(77.7),) = , g{m,m)=0, 
g{l,n) = l , g(777,m) = -1 , g{l,m) = , g{n,m) = 0. 

Remark 1.19 • The existence of a Newman- Penrose tetrad is insured by the existence of a 
global section of the fibre bundle of orthonormal frames: if (ej) (a = 0,1,2,3) is such a 
section, the following family of vectors: 



(1.2) 



is a normalized Newman-Penrose tetrad. It is obvious that a given normalized Newman- 
Penrose tetrad gives rise to an orthonormal basis of Tf2 with the following reverse fomulae: 

eg = -i=(r+n») = -ij(7n» + TO") 

e? = j^Q'-n'^) eg = ^(m»-m'^) 

• Because the structure of null cones will be considered later, we assume that a Newman-Penrose 
tetrad {l,n,m,rri) is given first and, in a second time, gives rise to an orthonormal basis (e^) 
(a = 0,1, 2, 3). 

• Up to an overall sign, there exist two unique spinor fields in £(§"^), denoted by o"^ and 
such that: 

I"- = 0^0^' ,n°- = L^T^ and m" = o^l^' . 
These two spinors are chosen such that the following normalization is satisfied: 

• There exists an alternative notation for this spin basis, which is consistent with the duality 
property used to describe spinors. We note, in S"* ; 

= and ef = l^. 
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We also introduce their dual spinors in Sa (e^,e\) which satisfy: 

e^£^ — , £a^o — 0; 



they are: 

e 



" — — and e\ = oa- 



The vector e° can be written in function of the metric as for a = 0, . . . , 3. The components 
of its spinor form g^^ , called the Infeld-van der Waerden, defined as: 



AA' _ a A A' 
9a — ^a^A^A' ■ 

are the coefficients of the decomposition of e'^ in the basis 

— 9a ^A^A' ■ 

It is then known ([I^, section 3.1) that the CUfFord muhiphcation of a Dirac spinor by the basis 
vectors can be written: 

Lemma 1.20 The Clifford multiplication of a Dirac spinor 4>a + by the vector e° is given by: 

el ■ {cl>A ® V^') = tV2g\A'^^' ® -«\/2.g=^^^>A 

Remark 1.21 ; The Clifford multiplication can be interpreted as a contraction with the correspond- 
ing vector of the basis (up to a factor ±i^/2) by writing: 



zl- {cj,A+i^^') - iV2g''''gi,AA'i^^' -iV2g--g^;- Pa 



ab AA' 

9h ^ 

V2g{ea, ea)gaAA''4'^ - iV2giea, ea)fff ^ 



As a consequence, the Clifford multiplication by the vector l^^ is the contraction with and 
conversely the Clifford multiplication by is the contraction by l^^ (up to a factor ±i\/2): 

1-{^A+ ^^') = tV2{nAA'^^' - . o^ 

n-i<PA+^^') = iV^ilAA'i^^' -l^^'M ^ ' 

We conclude this section by giving the abstract index expression of the Dirac operator on 4- 
spinors ([Hj, section 3.1): 

Lemma 1.22 The Dirac operator is decomposed as follows: 

bA + ^^') = iV2{Waa'^^' - V^^>a) 



1.2 Analytic requirements 
1.2.1 Distributions on spinors 

The purpose is to write weak sohitions for the Dirac equation. The theory of distributions must thus 
be adapted to ensure properties of symmetry for the Dirac operator and the Chfford multiphcation 
so that the construction of Friedlander can be used with few adaptations. 
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Fundamental properties The basic elements needed in the next section are sketched here. 
Spinor-valued distributions are defined in [3] to construct fundamental solutions for the Dirac equa- 
tion. They were also developed in (55] to construct a Fourier integral operator for the propagator 
of the Dirac equation. 

Definition 1.23 A distribution u on the set ©(Szjirac) of smooth Dirac spinor fields with compact 
support on fi, endowed with its usual Frechet topology, is a <C-linear continuous mapping from 
l^i^Dirac) to C, I.e. a mapping which satisfies for all compact K in £7, there exists a positive 
constant C and an integer m depending only on K such that: 

V0eD(S),|«(0)| <C||0|U,™,x 

The set of distributions on M will be denoted by T)' {^Dirac) and the duality bracket by <, >. 

Definition 1.24 The support of a distribution u is the complement of the largest open subset O of 
such that any smooth function (f> with support in O satisfies: 

< u,4' >= 0. 

The set of compactly supported distributions is denoted £'(SDirac) and is the topological dual of 

Dirac), SCt of SniOOth SCCtiouS of S Dirac OU CI. 

If u is a locally integrable section of S^j^,^^ = S"* ® Sa' , which can be written u — + r]A' , it 
defines a distribution by: 

Jn 

where the smooth section $ is split as: ^ — cpA + i'^ ■ 

We define now the action of the covariant derivative in a direction V and of the Dirac operator 
on distributions by: 

Proposition 1.25 Let u be an element of T)' {Suirac) cmd V be a smooth section nowhere vanishing 
of TO.. The distributions Vyu and pu are defined by: 

Vcj) e T)iSDirac), < VyU, (/) >b'(S*^^_^J,D(Sb,„,) = - < U^V V (p >D' {S*^^^J ,D{Sn,rac) 
V(/) e D(§Dirac), < Pu, >B'(S*,^^J,D(SD,rac) = ^ < P0 >I)'(S*,„J,D(Sb.„,) 

These definitions agree with the Leibniz rule and the fact that the connexion is compatible with 
the symplectic product on spinors. 

We also need to define the Clifford multiplication with a vector: 

Proposition 1.26 Let u be an element of Ti' {Soirac) o,nd V a smooth section of TCI. We define 
the distribution V ■ u in T>' Dirac) by: 

V</> e 2)(§£)irac), < V • M, >D'{S*^^^J.V{Smrac) = < U,V ■ (f) > D > {S o.rac) ■ 

Proof 

The representation of the Clifford multiplication is the same for the dual SJjirac- Consequently, 
if u = (/)^ +XA is in Soirac and v ^ pA' + 0"^ is in S^^i^ac tl^^n: 

< w, Ca • u >Si^^„^,So.™c= -iV^g'^AA'X^' - i^/^g"^^^' (t>APA'- 
We notice that this expression is symmetric in A and A' so that we can conclude: 
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Remark 1.27 When a distribution on Soirac is represented by a function from into Soirac, the 
symplectic product (•,•) on Soirac is used to apply the distribution on a section of Soirac- The 
duality bracket will be in that case written (•, ■)D'(Soirac),'D(Soirac)- 

The previous results need to be checked since the definitions given in (|1.25p do not work when the 
symplectic product (or the e spinor) is used. We first need the following lemmata on the action of 
Clifford multiplication and the Dirac operator: 

Lemma 1.28 For any $ and ^' Dirac spinor fields on Q and V a vector field on Q, we have: 

Proof : ft is sufficient to verify the result for an element Ca of the frame. We calculate (ca • 
in components. 

(ea = -tV2eA'B'9''^^'U0''' + iV2e^'' g'^AA'X^' Pb 

with <f> = + X^' and 'ii = pA + 9^' ■ Noticing that: 

AB a aB ^ aBB' _ aBB' 

£ 9 AA' -.9 A' = -£B'A'g ^ £A'B'9 

we obtain: 

(ea = -iV2e^''g''BB'Ue''' + iV^SA'B'g'"''''' X^' Pb = -($, ea • *).® 

In order to verify the symmetry of the Dirac operator for the symplectic product, we will establish 
the following lemma: 

Lemma 1.29 Let <f> and 5* two spinor fields on 17. Then we have: 

(p$,*) = ($,p*) - div{V). 

where V is a complex vector field. 

Proof : The formula is proved at each point of let then p be a point in fi. Let {fi) be a 
orthonormal basis on Q. such that, for alH in {0, 1, 2, 3}: 

V/, /i = at p. 

For this choice of basis, we have, at the point p: 

(P$,vl/) = (/,.V/.$,vI/) 

ie{0,l,2,3} 

= - E (V/.$,/.-*) 

i6{0,l,2,3} 

~ - ^ {V/;(<l>,/i • ^i) — {^,Vf^{fi ■ *))} (,) being compatible with the connection. 

ie{0,l,2,3} 

= - E {^/.(*' ■ *) - ■ V/.*)} since at p V f j, - 

ie{0,l,2,3} 

= (<i>,P*)_ ^ V/, ($,/,;. VI,). 

16(04,2,3} 

Introducing the complex vector field v defined, at p, by: 

3 



i=0 
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with fi —< fi,fi >, we notice that 

ie{0,l,2,3} 

is the divergence of V. 

We present an alternative way to perform this calculation with abstract indices; the Dirac spinors 
$ and \1/ are split on '^oirac'- 

$ = ® p^' 
We now lead the calculation in the usual way: 

4=0$,*) = e^''{lVAA'P^')4'B+eA'B'{~l){V^^'(^A)x''' 

= lyAA'ie^^'p^'i^B) + (-*)V-4^'(eA'S'</'AX^') 
+p■4'^Vf,^B-</'AH)V^,(x^') 

= iyAA'{p^'V) + l^''^'{c^AXA') 

-eA'S'p-^'«V^^VB + ie^^c^AyBB'x''' 

= ^VAA'{p^'i>'') + ^'^^^'{c|>AXA') + ^('J'-P*)- 

It must be noticed that, in this new calculation, the remaining term can obviously be identified as 
a divergence. @ 

Remark 1.30 The vector field 

3 

y = 5]f,($,/,;.xp)e. (1.4) 

1=0 

is encountered several times in the following. Though it is used to perform the calculation, it does 
not seem to he intrinsic. It is nonetheless easy to give a more intrinsic sense to this computation. 
Let us consider the complex 1-form uj on Q: 

TfK^C — > C 
h I — > {^,h-^) 

The dual vector of this 1-form is the vector p.4p . The calculation can then be easily reinterpreted 
when noticing: 

d*cu=l V/,($,/,-vI/)L, 

\4e{0,l,2,3} / 

★ being the Hodge dual and p the volume form associated with the metric g. 
Definition 1.31 Let u be in D' {Soirac), and X in C°°TVL. The applications defined by 

•irac 

and 

are distributions, denoted respectively by X ■ u and pu. 
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Proof : This is a straightforward consequence of the previous lemma and the Stokes theorem.!^ 

Remark 1.32 ; These definitions agree with the previous lemmata when u is in ©(Sui^oc)- 

From this point, aU the distributions wiU be assumed to be represented via the symplectic 
product. 

If / is in 'D'(M) and [/ is a smooth spinor field on il, we define the distribution fU by: 

V(/> e 1){§Dirac), ifU, <P)'D'{Sowaa).'D{SD,^aa) /' (U, (t>) >D'(K),D(R) ■ 

Composition of a function with a distribution In the following, the construction of distri- 
butions with support on a light cone will be required. One way to achieve this is to adapt the 
contruction of Friedlander in [5] in the case of spinor valued distribution. 

Definition 1.33 Let S he a smooth function on ^l, with non vanishing gradient on fi. 
Let f he a distribution with compact support on M. 
Then, the application 

(p e D(r!) ^ (fit), f cbip)VSip)Mp)] 
V Js{p)=t J 

where VS'(p)j/i(p) is the contraction of the measure on M with the gradient \/S (or the Leray 
measure on the hypersurface Sip) — t), defines a real distribution denoted fiS). This distribution 
coincides with the composition of functions when f is represented by a function. 

We need to apply this definition to calculate the action of the Dirac operator to a distribution on 

^Dirac 

of the form fiS)U: 

Proposition 1.34 Let f be an element o/£'(M), S a smooth function chosen as in definition \1.3S\ 
and U a smooth spinor field on M . 
Then, in the sense of distributions, 

p ifiS)U) = /'(5)V(5) • U + fiS)pU, 

where "S/ S is the raised gradient, i.e. 

V(5) :=^(Ve,u) e,. 

i 

Proof : Let $ g T^i^Dirac) and (/i) an orthonormal frame. is chosen with support in domain Vl 
where V /. fi are all zero. We calculate (p ifiS)U) , <I')D'(Soi,^^).D(SDi,ac) using the previous definitions 
and lemma [T. 291 

(P ifiS)U) , <i>)l,'(S„_),D(So_) = {f{S)U, P<i>)D-(So.™).D(So™) 

= </, / (C^,P*)MSt >£'(R),£(R) 

JS(p)=t 

= </, / (pC/,$)AiSt >£'(R),£(R) (1-5) 

JS(p)=t 

+ </, / V/. ([/,/;• $)/l5t >£'(R),£(R) (1-6) 

JS(p)=t 

where ^St is the Leray measure VS'j/i on the hypersurface St — {Sip) — t}. We calculate the two 
terms independently; by definition, (jl.Sp is: 

< /, / (pC/,$)M5* >£'(R),£(R)= (/(S')pf/,$)l,'(s„._),D(S^,_). 

Js{p)=t 
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and (|1.6p is calculated using the same idea as in lemma (|1.29p : 



< / ^niUJi ■ <^)flSt >£'(R),£(R) = < /, / div(u)^St >£'(R),£( 

JS(p)=t J S(p)=t 



where v is the vector field on O defined by: 

3 



V = 

i=0 



with fi =< fi,fi >. Noticing that: 

d 



div(u)^ = / div(w)/^Sj 

S{p)<t JSt 



and using the Stokes theorem 



div(u)/^ = / <\IS{p),v > 

S{p)<t J St 



JSt 
J St 



we obtain, accordingly with definition 11.311 

<f, (v^(c/,/^ • '^>)^ist >£'(r).£(r)- (/'(5)vs'- c/,$)i„(s^,_),i,(s„_) 



St 

so that: 

p {f{S)U) - f'iS)V{S) ■ U + f{S)pU ® 

Spinors and bidistributions Keeping in sight that the purpose is to write an integral formula 
(or representation formula) for a Cauchy problem, we must be able to apply twice a distribution to 
spinor fields. This is what bidistributions are made for. 

We define the product Kl of two smooth sections of Soirac with compact support by: 

25(SDirac) X 2?(S£)irac) > '^{^Dirac) ^ "^{^Dirac) 

("I>,^') I — y e r2 X 17 ^'(p) (g) $((7)) 

where (S" must be understood as tensor product of spinors in different variables. The vector space 
generated by these products is denoted by T){§uirac) ^ '^{^Dirac)- 

Definition 1.35 Let u and v be two distributions in D' {SDirac) ■ The bidistribution uMv is an 
application from ^{SDirac) ^ "^{^Dirac) defined by, for every (0, ■0) G ^{SDirac) x "^(^Dirac): 

The vector space D' {S>Dirac) ^ 'D' (SDirac) generated by these products is called the space of spinor- 
valued bidistributions on Soirac- 
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If (j) is in D[E>Birac) and m is a spinor bidistribution, then is still in D' {E>Birac)- It can conse- 
quently be still applied to a function in ©(S/jirac)- 

A special type of spinor valued distribution that will be encountered in the following is the Dirac 
distribution. 

Definition 1.36 We define the Dirac distribution (or Dirac mass) in p' , denoted by 5pi by: 

e 2)(§), (V, = </>(p')- 

It must be noted that this distribution can be written in the form t{j)' ,p)5pi ([6], chapter 6) where r 
is a linear transformation from D{EiDirac) in the variable p to D^^^irac) in the variable q satisfying 
'''{p' tP') — koiraa ^nd cau consequently be written as: 

V(/) e Ti{S), (t(p',_p)V, 0)p,{I)'(Sd„.„c).d(Sd.™<:)} = 
the duality bracket being computed in the variable p. Since 

£ ^A^B — ^ ana Ea'S'^q' — ^-L 

and t{p,p) satisfies: 

ij{p,p),(j)(jp)) (j){p) 
it can be explicitly calculated at p = p': 

t{p,p) = -e^' Me^' +£^' Me^' +elme\-e\Mel 

= -t'^' mo^' + 0^' mi^' - OB^LA + iB^OA. (1.7) 

Such a function t is chosen explicitly later (see equation p.9p ). 
1.2.2 Fundamental solutions of the wave equation 

We now apply to the spinorial wave equation the analytical tools used by Friedlander in 6 for 
the tensor wave equation. An alternative method has been used by Klainerman and Rodnianski 
to construct an approximate fundamental solution in [14] . Though their method is more flexible 
and well-suited to obtain estimates, it is not appropriate here since, as we will see, the regular part 
(the tail of the fundamental solution) is needed to write down a fundamental solution. V. Moncrief 
used Friedlander's method in a paper with D. Eardley ( 4J) for the Yang-Mills equations in the 
Minkowski space and for the Maxwell wave equation in [16 on a curved space-time. 

We first consider the spinorial wave operator p . The Schrodinger - Lichnerowicz - B5chner 
formula gives that for any in D(§£)irac): 

p^0 = n(/)-f iscal(/i (1.8) 

where □ = — V^V-' . Since the index notations are used from the beginning, a index version of the 
formula with its proof is given: 

Proposition 1.37 (Schrodinger-Lichnerowicz formula in index version for spin i) 

Let (j)A be a smooth section of Sa ■ Then we have the following relation: 
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Proof : The reader should refer for intermediate results to [1T]{4.9.2 and 4.9.17). 

= ^e'^^VHH'y^^'eBC^A + ^Scalg^s (formula 4.9.17 in [H]) 

= ^VcC'V^^Vb + lsc'd\g(pB® 
Z o 

Remark 1.38 • This version agrees with the previous one when noticing that the operator 
VsA'V'^"^ is in fact, due to the renormalization induced by the Clifford multiplication, the 
projection on §b of l/2p^. 

• A generalization of this formula to arbitrary spin is given later in subsection \3.1[ 

Since is a geodesically convex domain, it is possible to define globally on Q the squared-distance 
function: 

where 7 : [0, t] ^ 51 is the unique geodesic from p to q. 

To write the fundamental solutions of the wave equation, it is necessary to construct distributions 
with support on a cone: using definition 11.331 let us consider the distributions 

S^iTpiq)) and H^{Tp{q)) 

where S is the Dirac mass and H the Heaviside function. These distributions have support respec- 
tively, for p fixed in fl, in C^{p) and 3^{p)- 

Remark 1.39 It is important to notice that these distributions do not satisfy definition \1.33\ since 
the gradient ofTp{q) vanishes at the vertex of the cone. Nonetheless, considering the distributions 

S^{rp{q)-e) andH^{Tp{q)-e) 

with e positive avoids the problem. The results can then be obtained using a limiting process when 
e tends to zero. This method will be used later to expand equation (|2.8p . 

It is known that the operator admits fundamental solutions ([H],[S]): 

Theorem 1.40 There exists two bidistributions on il, G^[p) that satisfy: 

y{p,q)^n\t>lG^{p)^5,{p) 

in the distribution sense. These two bidistributions can be written: 

G±(p) = U,{p)5^{T,{p)) + V,{p)H^{T,{p)). 

where U and are smooth functions of the variable {p, q). q being fixed in il, the support of G^{p) 
is then in Q^{q). 

The structure of the fundamental solution obtained by Friedlander is the following (the reader 
should refer to [5] for more details.) 

a. The function U in the singular part can be decomposed into two parts, Uq{p) — kq{p)fq{p) 
where: 
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(i) the bispinor fq{p) satisfies: 

VTq(p)Vif<,(p) = and fp{p) = Tp{p). (1.9) 

This equation can easily be reinterpreted as parallel transport in the variable q of the 
bispinor identity along the geodesic from p to q. 

(ii) the function kq{p) satisfies the transport equation: 

2 < Vr,(p), Vkgip) > +(nr,(p) - S)kg{p) = O and kp{p) = ^. (1.10) 

kg (p) measures the difference between the measure induced on 6+ (p) (1 C~ (g) and the 
measure on the standard sphere S"^ in the sense that, if p is in the future of q: 

/^e+(g)ne-(p) = kg{p)r^ns2 

where Me+(g)ne-((jj)) is the Riemannian volume form induced by the metric g on C+(g) fl 
C~ {p) and /i52 the standard volume form on the two dimensional sphere. 

b. The regular part of the fundamental solution can be obtained by solving the characteristic 
Cauchy problem: 

aVg{p) = OforpGa+(g) 
Vg{p) = (p) for pee+{q) 

where Vg{p) satisfies the transport equation: 

2 < vrg{p),vv°{p) > +(nr,(p) - 4)i/"(p) = -d^u. 

For later convenience, the fundamental solution must be split over the decomposition of the 
Dirac spinors: 

G^{p)='G^g{p)+'Gi{p) 

A B 1 

1 

The notation A means that the part of the bidistribution in the variable q acts on spinor fields in 
§A- Their fundamental part is denoted by, respectively, [p] and ^Ua'(p). 

A B 9 

Two backward and forward fundamental solutions for the wave equation can then be constructed. 
For Dirac spinors, these fundamental solutions are the distributions: 

D^G^{p) 

In terms of indices, these fundamental solutions are written: 

(p) on §A El and V%b'^Ga'{p) on §^ K § 



-q ' 
A B 



B- 



Finally, we state the following theorem concerning the existence and the structure of the funda- 
mental solution for the Dirac equation for Dirac spinors. 

Theorem 1.41 There exist two fundamental solutions for the Dirac operator p, G^{p), with sup- 
port in G^{q), for q fixed in fl, such that: 

y{p,q)Gn\p'G^{p)=6g{p) 

in the distribution sense. These two fundamental solutions are obtained by applying the Dirac 
operator to the two fundamental solutions of the wave equation: 

G^{p) = p'G^{p). 
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2 Derivation of the integral formula for Dirac spinors 



This section is devoted to the derivation of an integral formula for Dirac spinors for the characteristic 
Cauchy problem with data on a future null cone. In this context, we will work with the forward 
fundamental solution G^lp) which will be denoted with no ambiguity Gq{p). The singular and 
smooth parts of the forward fundamental solution for the wave equation will be denoted Uq{p) and 

The point po being fixed, let p he a point in the future of po in Q. We define, for these two 
points: 

. aip) = e+ipo)ne-{p) 

• D(p) = e+(po)nr(p) 

• §{p)^3+{po)ne-{p). 

• v(p) = a+(po)nr(p) 

Since is geodesically convex, these instersections are well-defined (in fact, the hypothesis of global 
hyperbolicity suffices). 



2.1 Representation formula 

The first step to obtain a representation formula is to solve the problem with source: 

pu = f. 

The following lemma is a direct transcription of lemma 5.5.1 in [5j: 

Lemma 2.1 Let f in £{SDirac)- 
Then the distributions defined by: 

V0 e D(§z5.™c), {u, 0)p := (/, (G±, 

are solutions of the problem: 

Proof : The calculation is made first formally. The justication of each step will be carried out later; 
it will be sufficient to check that each duality bracket is well-defined and that all the operations 
involved (symmetry on Dirac operator, . . . ) are legitimate. 
Let 4> be in ViSoirac)- 

{p^u,^)p = {u,p''cl,)p (2.1) 
= (/,(G±,pV)g)p by definition of u (2.2) 
= (/,(/)) by definition of G±. (2.3) 

It must be checked to insure that (12.21) exists that the function: 



P^(G±,p^ 

is smooth; we have: 

(2.4) 



(p^G^,]Z)V),-(G,^(r)^0),- / (C/p*(g),(p')^0)Mr,(,)(g)+ / {V^{ql{WYmq) 

Je+{p) Js+ip) 
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where fJ^Ppiq) is the Leray form associated with the function Tp{q), i.e: 

Mrp(g) = VTp(g)j^. 

Let TT : he a chart recovering ^l (which exists since ft is geodesically convex). The image 

by TT of p and q are respectively denoted by y and x. There exists a diffeomophism ^ — > a; = h{y, ^) 
from 7r(17) into W^, where ^ — (C*^, C^, C^, C'^) is a coordinate system centered at y, Minkowskian in 
q and such that the vector (1, 0, 0, 0) is timehke and future oriented. In this coordinate system, the 
measures /x and Mr,(p) are expressed as: 

m) = k(y,Odt. and ^r,(p) = k[y,^)- 



with d^ = dC° A d^i A d^2 ^ g^jjjj 



c+(g) = {eie" = Vier + ier + ier} and = {m" > vwrwrw}- 

The integral p.4p can then be rewritten: 

(2.5) 

(F4^)(y),((P^)20)(%,O)fc(2/,C)de (2.6) 



C">\/(«')'+«')'+(«')' 
which is clearly a smooth function of y = n(p). 

Since / is a distribution with compact support, there exists K' , an integer N and a positive constant 
C such that, for any smooth function ■0, the following estimate holds: 

|(/,0)|<c J2 \K^o^-\y)\\- 

Let A" be a compact of fl. Assume that cj) has its support in K . Then the previous inequality gives 
\{Li^)\<C'Y^ sup ||5^"(G.-.(,)(^-i(x)),0),||. 

Using the expression of (G^-i(j,)(7r^-'^(x)), 0)2,, its derivatives 9^(G^-i(y)(7r^^(a;)), (^)2; are bounded 
by the derivatives of on K: 

snp\\d;;{G^-i^y^{7T-\x)),<PU\<CK,K',a J2 ™P 

where the constant Ck,k' .a is determined only by the derivatives of U , V, h of order up to fc + 1 
on the compact K x K' and its image by tt. Finally, we obtain: 



which means that equation (12. 2p is well-defined. @ 

Let u be in £(§£)irac)- The following proposition gives a representation of u in term of its data 
on a null cone: 
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Theorem 2.2 Let u be a function with future bounded support. Let po in fl. Then, we have in the 
distributional sense: 

uH+iToj^P" (( f {p^u,Up{q))^^rM+ f {P'u,Vp{q))t,rM 



+ {VPTo-u,Up{q))firo,rM+ {V^To ■ u,Vpiq))fXroip) ] H+iVo) 

Ja(q) JD{q) ) ^ 

where the two-form /^ro.r, (p) *s obtained via the factorization: 



V0 e C^{n X n), / 0/iro(p) A Mq) = / / ^t^T„,rM a fliq), 
where A is the set {{p,q)\p G C+(po) o,nd q G G^{p)}. 

Proof Let w be a function with future bounded support, that is to say that the intersection of 
supp(u) with any future null cone is compact . We use here the property of the fundamental 
solution with lemma [^TT] with / p^{uH^(ro)): 

{uH+{To),cl,)p = {pP{uH+{To)) ,{Gp,<j>)^)^ (2.7) 

{uH+{To)) , (p^Gp, 



p^(ui/+(ro)),(Gp,pV) ) (2.8) 



The duality bracket (12. Sp is properly defined since the function p t-^ {Gp, p'^4>)q is a smooth function 
with support in the future of Supp((/)), that is to say Uqgsupp(0)3^('?)5 ^^id since u has future bounded 
support. 

The duality bracket (|2.8p is then developed. The first step consists in differentiating the dis- 
tributions MiJ"'"(ro). As already noticed in remark ri.39|, the distribution uiJ+(ro) is not of the 
type given in proposition 11.341 since VFo vanishes at p^. To avoid this difficulty, we consider the 
distributions uH'^{To — e), where e is a positive number. This derivation gives, since proposition 
11.341 can be applied: 

p {uH+iVo - e)) = ipu)H+{To - e) + VFo • uS+{To - e) 
which becomes, when e tends to zero: 

p {uH+{To)) = ipu)H+{To) + Wo • ud+iTo). 
The bracket (|2.8p is written as the sum of four integrals: 



fp^(«i7+(ro)),(Gp,pV) ) = / / (p^^,(?7p(<7),P»)A*r,(g)AA*(p) 



3+{po) Je+{p) 

(P^u, {Vp{q),p''^))fi{q)Afi{p) 

3 + (po) JS + ip) 

(Vfro • u, (C7p(g), pV))/^r,(g) A /.r„(p) 

e+(po) Je+{p) 

{VPTo-u, {Vp{q),p'^c^))^l{q)^^lrM. 

e+(po) J3+(p) 
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where /ipo and /ipp are the Leray measures associated with Tq and Tp respectively. 
Switching the order of integration of the variables, we get: 

(uiJ+(ro)),(Gp,pV) ) = / / ((p^u,t/p(g)),pV)Mr,(p)AM<z) 
V ^ J 3+ ipo) J e- {<])ns+ ipo) 



3 + (po) J3-(q)n3 + ipo) 

f (V^To • u, Up{q)), pV)A'ro A ^i{q) 

3+(po) Je-(q)ne+{po) 

[ ((VP ro • u, Vp{q)), pV)Mro (P) A fi{q). 

3+(po) "'3-(9)ne+(po) 



Finally, the duality bracket (|2.7p is: 



(«i/+(ro),0)= / (pVC/p(9))Mr,(p)+ / {p'u,Vp{q))firM 

\\J8{q) Jv{q) 

+ f {yPro-u,Upiq))firo,rM+ f (VPro-u,yp(g))/iro(p)W+(ro),PV 

Ja{q) Jv(q) / j 
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9. 



which means that, in the sense of distributions, u satisfies, using the symmetry of the operator p 
uH+{To) =P'(l f {p''u,Up(q))t,r,{p) + [ {P'u,Vp{q))firAp) 

\\J8iq) Jv(q) 

(V'To-«,f/p(<7))Mr„.r,(p)+ / {W^T^ ■ u,Vp{q))^lrM] H+{U)\ M 

cr(q) Jv{q) J J 

A direct application of the previous theorem is the first integral formula for the characterictic 
Cauchy problem: 

Proposition 2.3 Let u be a smooth solution of: 

pu^O 

Then u can he expressed in 3^(po) *^ function of its restriction to the cone C^(po) ^U' 

"la+(po) - PM ( / (VPro-ii,c7p(g))Mro,r,(p)+ / {V^V^ ■ u,Vp{q))^irM] H+{tM . 

\\Ja(q) J-D{q) J J 

Remark 2.4 a. This formula is not the final stage of our calculation; the fact that it only 
depends on initial conditions will he stated later. This is the purpose of the next subsection. 

h. The vector VFo being null along the cone &^{jpq), Clifford multiplying with VFo means in fact 
contracting with the spinor form o/ VFq; a direct consequence of this is the fact the Clifford 
product of the 4^- components Dirac spinor s with VFq does only involve the two components, 
Mo and . The two remaining components are recovered using the constraints equations ( cf. 
lemma \37?\ below). 



2.2 Integral formula 

The integral formula is derived in three steps: 
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• construction of the appropriate geometric tools (derivation of measures, spin basis); 

• interversion of the integral and the Dirac operator; 

• and finally obtention of an expression of the singular part in terms of geometric quantities 
and initial data. 

2.2.1 Geometric data on the cone 

This section is devoted to the calculation of the relevent geometric quantities for the intersection 
of C+(po) n G^{q) = cr{q) for a given point q in the future of po- This is widely inspired by section 
4.14 of [21] ■ There are also some calculations of interest in the work of Frittelli, Newman and al 
([7], for instance) and Nurowski - Robinson([T5]). This kind of calculation is also very common in 
the study of Ricci flows. 

We first choose a parallely transported vector field / along the null cone C+(po): 

Vil = 0. 

We then consider, for a given point q in 3^{po), a point p in a{q). We construct at p a Newman- 
Penrose tetrad: 

a. the first null vector is the vector l{p) at p; 

b. n{p) is chosen on the future oriented null geodesic from p to q such that g{l, n) — 1; 

c. we complete the basis by taking a pair of complex null vectors m{p) and rn{p) in the orthogonal 
of the vector space generated by (/, n) such that g(m,rn) = —1. 

A Newman- Penrose tetrad is then obtained at each point q' on the cone C~ (q) : let p' be the point 
in a(q) lying on the unique null geodesic from q' to q; the Newman-Penrose tetrad is obtained in q' 
by parallely transporting the one at p' along the unique null geodesic from p' to q' . 

Remark 2.5 This construction cannot be realized globally on the intersection {po)r\C^ (q) = cr{q) 
which has the topology of §^ . It will be necessary to make this construction on two different open 
sets and then glue these constructions together to obtain the result which only depends on I and n. 
We assume then that the construction is done on one open set. 

This choice of Newman Penrose tetrad gives us: 

• a basis of Til C and, consequently, up to a sign, a spin basis of §>-^ that will be denoted by 
(0^,6^); 

• if g is fixed first and p is chosen on a{q), the vectors m and m span the tangent plane to a^q) 
at p: 

TpCr{q) = {Am + Am|A G C}; 

due to obvious topological obstructions (see remark [^5]) . this construction cannot be extended 
globally to all a{q). 

• the choice of Z, which is parallely transported along the generators of C+(po), and n, which 
is parallely transported along the generators of 6^ {q) , gives rise to two affine parameters tq 
and r along the null geodesies on these two cones. 

• these two affine parameters give rise to two parametrizations by the sphere S*^ of a{q) using 
the exponential map at po and p respectively: 

exppo : — > n 

oj I — > exp (ro(w)w) 
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and 

expp : 

uj I — > exp (r(a;)ijj) 



expp : — > ri 



Let g be a point fixed in 3^(po)- We consider a point p on cr(q). In a neighborhood of p, on 
cr((3'), is defined a Newman-Penrose tetrad {l,n,m,rn). The dual basis in Tfi* (g) C is denoted by 
(L, TV, M, Af ) for which tlic following lemmata are true: 

Lemma 2.6 The induced metric on <t{p) is ~2MM, the volume form A M and the mean 

curvature vector: 

H ^ 2{p'l + pn) 
where p and p' are the real spin coefficients: 

p = —{I, Vmm) and p' — -~{n, Vm^) 

Proof : These results are straightforward consequences of the presentation concerning two-surfaces 
in (section 4.14, proposition 4.14.2 sqq.) 

The reality of the spin coefficients is stated in proposition (4.14.2) of [H], whenever I and n are 

orthogonal to a spacelike 2-surface (here cr{q)). 

Since (m,rn) span Ta(q), the second fundamental form is: 

V(X, Y) e Ta{q),II{X, Y) = {VxY, n)l + (V^F, l)n 

so that the mean curvature vector is: 

H = —(II(rn,m)+II(rri,rn)) 

= - ((((5'm, n) + (6m, n))l + {{5'm, I) + {5m, l))n) 

Since (see [H] (4.5.28) together with (4.5.29)): 

6'm — (/3 — oi)m — 'p'l — pn 
Srn = (a — (3)m — p'l — 'pn 

and since p and p' are real, we obtain: 

H = 2{p'l + pn).% 

In order to compute the Leray forms associated with the distance function, we use the expressions 
of the gradients of the distance functions Fq and Tq-. 

Vfro(p) - 2rol{p) and V^T^ip) = 2rn{p) (2.9) 

The Leray forms can then be expressed using the dual basis of the chosen Newman-Penrose basis. 

Proposition 2.7 The Leray forms p-Vg, /^r, and pvoX 

pLTo = — ^ AM AM 
2vro 

Air, = —LAM AM 

1 — 1 

itr^r Ar^r ^ ' 

where VPro(p) = 2rQl{p) and VPTq{p) = 2rn{p) 
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Proof : The volume form on O can be expressed in terms of the Newman-Penrose tetrad as: 

H = -LAN AM aH 

so that, since dFo = 2roL and dF, = 2rN, we obtain immediately: 

ur„ = ^J—N AM AM and ur = -^L AM AM. 
" 2iro ^ " 2ir 

The calculation of /Uro,r is obtained through the factorization given by Fubini's theorem: 

V0 G C^in X ft), [ 0/xro(p) A Mq) = [ [ '/>m,r,b) A fi{q), 

J A J3+{po) J<y{q) 

where A is the set {{p-,q)\p € &'^{po) and q€:&'^{p)}. We get the desired expression of /iro.r,: 

Mro.r, = — ^MAM= 

War 4ror ^ ' 

The next step consists in determining the variation of the metric when q is in J+(po)- We 
first establish the technical lemma: 

Lemma 2.8 Let (J^,h) be a smooth semi-riemanian manifold with metric h and Levi-Cevita con- 
nexion D; let X be a smooth vector field on 3Nf. Let (Mp,^) be a submanifold ofJ{ such that the g 
metric induced by h is non degenerate and depending smoothly on a parameter p in'N in the sense 
that there exists a smooth manifold S and a smooth map / : K x S — > K which satisfies: f{p, -k) is 
an immersion and f{p, E) = JAp. 
We denote by the induced volume form on Mp. 
Then: 

I^xl^p = -h{H,X)^p 
where H is the mean curvature vector field on Mp. 

Proof : The Levi-Cevita connection induced by g on Mp is denoted V. 
Let p be a point in K and q a point in Mp. We choose around q a map 

(V, ix\x'^, . . . ,a;"), (a;"+\ . . . ,a;"+'=)) normal at g and Mp n F = {a;"+i = • • • = a;"+'= = 0} such 
that, at q: 

Va^.a,. =0. (2.10) 
The volume form on Mp around q can be expressed: 

/Up = \det{gij)\^dx^ A dx^ A • • • A dx". 

We calculate the derivative: 

Dx^ip = Sign(det(ffi,))g^^|det(5i,)|^da;i A da;^ A • • • A da;" 
and then evaluate at q, where the coordinate system is normal: 

n ^ 

^xfip = -£iDxg{d^i,d^^)iXp 

i=l 
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with Ci — g{d^i^ d^i). Since the connection D on is metric, we get: 

n 

i=l 
n 

n 

1=1 

Since Dq .^^.i — Vg .d,j.i + II{dxi, d^i), we finally obtain: 

n 

We then notice that, for i in {tl, n}: 

h{[X,^x^],^x^) = -£i^r- and Dq .h{X,d^^) = 



ox^ ^ ax* 

Since Va . d^i is vanishing at q, the only remaining term is: 

VxA*p = -h{X,H)^ip.<^ 

A straightforward application of this lemma is the proposition: 

Proposition 2.9 Let f : fl^ ^ ^Dirac be a smooth mapping. 
Then the following formula holds: 

f f{q,p)t^.i,){p) = I P'f{q,p)+S/Pro-i7ff{q,p)-2p\/'^ro- f{q,p)fi,^,){p) 

Jcr(q) Jaiq) 

where ro, being a function of both p and q, satisfies VFo — ^r^l. 

Proof : Let V a vector field on O. We work with the exponential map centered at po- cr(9) can 
then be parametrized by S"^: 

w 1-^ exppJro(g,w)w). 
Let us consider the variation of a{q) defined by, for some positive e: 

]-e,e[xS'2 ^ f7 

{t,uj) ^ expp^(ro(g + iV",tj)w) 



Since 



V«.(/(g,p)) = V^(/(g,expp^(ro(g,a;)c^)) 



and 



V?./(g,p)+(^/(g 



this gives, using lemma 



V'v^ / /(g,pK(,)(p)- / V?,/(q,p)+V?,roVf/(g,p)-2pVt.ro/(g,pK(,) 
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so that, when choosmg an orthonormal basis (e^) on Q, we obtain: 



Jij{q) Jryiq) 
We finally establish the following proposition: 



Proposition 2.10 Let f : il^ ^ ^Dirac be a smooth mapping. 
Then the following formula holds: 



Pi f{q,pWo{p)^ PV(9,pKo(p)+ / V%-/(g,p) 



, Ai<T(g) {P) 



2r 



where tq and r, being functions of both p and q, satisfy V^Fq — 2ro/ and V^F^ = 2rn. 

Proof : We use exactly the same method as in the proof of proposition [231 Using the paramctriza- 
tion of the exponential map centered at pq, we have: 

rro(<3,w) 



P / /(g,p)A*ro(p) = P fiq,rio) — k{u,r)dus2 

JD(q) JS'^ Jo ^'^ 

= / Pfiq,pWM+f V%./(,,m'^"''^'"^^ 



Pf{q,rLo)iirM+ ^'ro-f{q,p) 



2r 

l^'y{q)iP) ^ 



Dig) J<T{g) 2r 



2.2.2 Derivation of the integral formula 

We now consider the characteristic Cauchy problem on Q: 

pu = Oon3+(po) 



u = 9one+ipo) ^^'^^^ 

where is a smooth spinor field on C^(po)j whose support does not encounter the vertex of the 
cone and satisfies the constraint equations given by lemma 1X71 

Remark 2.11 The term "spinor field on the cone" must be understood as "trace on the cone" of a 
Dirac spinor field on i7 and not as a spinor field constructed as spinors on the manifold C~'"(po)- 

The basis constructed in the previous section is used to split the spinors: 

e = e'e^' + ci£^ = e'o^' + e'-^^' + co(-m) + cio^. (2.12) 

u will be split on S^i (g) S'^ : 

U ^ (t)A + 



The solution of (|2.1ip can be written in function of its data on the cone and the basis {o^, l^)- 



27 



Theorem 2.12 Let u be a solution of (|2.1ip . Then, for any q in 3^{po)' 

•I (y{q) \ ' / 



J<y(q) 

ipPV„VPTo■u)^lr„{p) 



Remark 2.13 It is possible to obtain a representation formula for the Goursat problem for the 
Weyl equation: 

V^^'(t>A = 

by projecting the solution obtained in theorem \2.12\ on the subspace of Dirac spinors Sa- 

Proof : Let g be a point in 8^{po)- Using proposition 12. 3| proposition 12.71 and proposition 12. 10[ 
we have: 



(Vfro-^/,p«Fg)/xr„(p). 

D(g) 

The bracket in the first integral can be calculated as follows: 
(^^' ^2rro " kp{q)(^Tp{q),^n-u 



iV2kp{q) (^Tp{q), (^-0^0^ (pA + o^oa'V' 
i^/2kp{q) (^Tp(g),^ (-000-^' +i/;1'oa)^ 

iV2kp{q)j- (^~Mp)o^' (q) + i^^' {p)0A{q) 



28 



We then use proposition 12 . 91 to calculate the first integral: 



zV2p' n kp{q) (-Mp)o^\<i) + i^^' {p)oA{q)) 



Ja(q) \ 2r 2r J 

-ipVVo • (^-Mp)o^' {q) + {p)oA{q)) f^a{q){p)- 

In order to simplify the calculation, the previous formula is projected on Sa- The singular part on 
this subspace is written, after expansion: 



A = iV2 f iV2V%, I kp{q) i 

Ja(q) \ \ 



^Mp)o^'i<i) 



2r 

■A' I 



Mp)o (g) 

2r 



-{^V2)k,iq)pV%,ro ^^^^^ (p) 



(9) 



r 



-2kp{q)pV\^,ro I 1 Pa{q)(P)- 



Mp)^\a> ( kM { + v^^.ovr ( kM { M^)^' 



r 



~2kMP^\A>r. f '^"(^f"'('^) 



Expanding all the products: 



A 



{MP) {^o^\q)^\A> (^) + (^) ^\Aro-\q)) 



r 



V^cf>o + V^['^]<j>o]iV\^,ro)o^iq) 



2 (^tM^ p(Vi^,ro)0o(p)o-''(g)) M.(,)(P) 

1 /ytPj, /'^\ fY7l „ \^A' 



f 



Mp)yAA'(^i^'^y'''iq?jfMq)iP)- 
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Since the quantities which appear in the integral are the restriction of u and its tangential derivative 
along the null cone C"'"(po)j 4> can be replaced in the integral by the data of the Goursat problem 
Co: 

Ja(q) \ r / 

^/(,) ((^) "■''('^O ^-(^)(^)- 

We obtain the complete formula for Dirac spinors by adding the corresponding quantity on S'^ , 
meaning: 

J (j(q) \ ^ J 

Noticing that the calculation has been done for (C/p(g),Vro • u) in order to use the definition of 
t{p, q), we obtain the complete formula using the antisymmetry of the symplectic product. @ 
It is now possible to obtain the formula established by Penrose in [3D] in the Minkowski case: 

Theorem 2.14 (Penrose) Let u be a solution of (12. lip . 
Then, for all q in 3^{pa), u can be written: 

<q) = ^^^^ - 2P^''(^')) (^r'^o) OA(g)/i.(,)b) 

+ / i(^fCo(p) - 2pCo(p)) {V\^,r,)o^' {q)p.i,){p) 

Remark 2.15 First of all, the meaning in the context of a flat space of the choice of the basis 
constructed in the previous section should be made precise: 

• the spinor is chosen to be constant on the null generators of the cone; the affine parameter 
rQ is measured with respect to the vector — 0^0"^ ; 

• a direction on the cone C^(po) being given together with a point q in £j^(po); ^ef p be the 
intersection of C^((?) with this direction on the null cone from po; the spinor is chosen so 
that n" — L^T^ is colinear to the vector pq and satisfies: oaI'^ — I7 the affine parameter r is 
measured with respect to the vector n" — l^T^ ; 

• the basis is completed by the two vectors — o^T^ and m" = l^o"^ . 
This construction is the "fiat" version of the one made using parallel transport. 
Proof : As done in 20J, it is sufficient to remark, for a direction uj on the cone C+(po): 

q = Po + rol^iuj) + rn°-{q, u), 

which implies: 

VV = r, V% = n" 

and kp{q) = 
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Remark 2.16 It is interesting to note that the term that carries the curvature information in the 
singular part is: 




(2.13) 



It is somehow difficult to give a precise geometric interpretation to equation ()2.13p . Nevertheless, 
clues can be found in theorem 4.2.2 in '61 that states that {k/r) measures the growth rate of the 
measure fJ.a{q) ■ 



3 Generalization to higher spin 

In this section, we obtain an integral formula for solutions of the Goursat problem for the Dirac 
equation with arbitrary spin. The derivation of the formula is based on the representation formula 
for the Weyl equation which can be extracted from theorem [2321 

Let us consider the characteristic Cauchy problem for spin ^ — s > 1 {n being the number of 
indices of a spinor): 

( \7^'^'uAB...F = Oong+(po) -^x 
1 U00...0 = 6'oo...o on e+(]3o) 

where uab...f satisfies the symmetry conditions: 



UAB...F = U(^AB...F)- 

First of all, it must be noted that, on an arbitrary curved space, the problem p.ip cannot be 
set if a consistency condition on the conformal curvature is not satisfied ([T], [S] and [I7j for the 
Rarita-Schwinger case for a treatment of the Cauchy problem). It is known that for the Dirac 
massless equation for low spin (n < 1, i.e. scalar wave, Dirac- Weyl and Maxwell equations) this 
condition is always satisfied. For higher spin, it is satisfied whenever the space-time is conformally 
flat. Nonetheless, it is expected that the method could be adapted to the Rarita-Schwinger case 
which requires the space-time to be Ricci flat. 



3.1 Generalization of Dirac equation to higher spin. 

The construction that was made before for Dirac spinors is adapted here to spinors of higher valence 
so that the symmetry conditions of the Clifford multiplication and Dirac operator still hold. 
Let us consider E the fibre bundle defined by: 

^ = Sab...f®S^'g...i- 
This fibre bundle is equipped with the symplectic product obtained from e: 

and a Clifford multiplication by vectors: if u — (I>ab...f + V'"^ q j belongs to E, we define Ca • u, 
where (ea)a=o,...,3 is the basis constructed in subsection ll.il 

Ca • M = ~iV2g'^^^'(PAB...F + « ^2^'^^^, V''^ g. . ./ ■ 

We finally define on smooth sections u = <j)AB...F + g of E the following operator (that 
will be denoted by p as the Dirac operator for Dirac spinors): 

pu = iV2{ - V^^'(j)AB...F + Vaa'V-"^ G.../)- (3-2) 

The distributions on smooth sections of E are defined using the (non degenerate) symplec- 
tic product £ in the same way as in section 11.21 The duality bracket will still be denoted by 
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(: )i3'(E),D(E))- Let u — (pAB...F + i'^ q j and v — £,ab...f + C"^ g / be two smooth sections of E. 
We have: 

We finally extend the Schrddinger-Lichnerowicz formula to arbitrary spin: 

Proposition 3.1 (Schrodinger-Lichnerowicz formula for arbitrary spin) 

Let he ipF...i o, smooth section ofE>F...i (n indices). 
Then the following formula holds: 

VbA'^^^^'^F...! - ^Vhh'V^^VbG.../ 

^ F^'4'DG...I — Q^->jjFD...I — ■ ■ ■ — Xg^ j^->pFG...D 

where Xabcd is the curvature spinor: 

XaBCD = -^RaX'B^ CY'D^ , 

R = Rabcd being the Riemann curvature tensor of CI. 

Remark 3.2 It must be noted that the potential of the operator p^, though linear, is no longer 
scalar and not even symmetric. 

Proof : the proof is almost the same as the proof of proposition 11.371 

VbA'V^^Vf.../ = e^^VBA'^c^F...! 

= (V[s|A'Vj^]^F.../ + V(s|A'Vjt^)V^F.../) 

= le^'^^NH'^^^'sBC^F...! + e^^V[B\A'^fc]^F...I 

= ^Vhh' V^^VsG.../ + e-^'^V[B|A'Vj^]V'F...7 

The spinor ipF...i is then split as the sum of tensor products of spinors of valence i, and, then as 
explained in [H] (vol. 1 p. 245, together with formula (4.9.4), (4.9.5) and (4.9.8)), using the fact: 

V[_B|A' V|^jU_D = —Xg(j^ DUE 

for any smooth section of (formula (4.9.8) in [21]), we obtain: 

V[s|A'Vj^]'(/'F.../ ~XgQp^^DG...i — Xg^Q^^FD...! - X iI:fg...d 

and finally: 

Vsa'V^^Vf.../ = -\vHH'V""'i^F...I 

- e^'-^ [Xg(jp^ipDG...i - Xg^Q^^pFD...! - Xb(jj^i};fg...d) 

= -\vHH'V""'i^F...I 

~ Xg^ p^1pDG...I — Xg^ Q^ipFD...! — • • • — Xg^ j^^FG...D® 

As an obvious consequence of the definitions chosen for the Clifford multiplication and the Dirac 
operator on E, the following proposition holds: 
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Proposition 3.3 The Dirac operator p on E and the Clifford multiplication by a vector field v on 
fl are respectively symmetric and skew symmetric with respect to the duality bracket (, )d'(e).d(e) 
that is to say, for any 4> and ip smooth sections of E with compact support : 

(0,P'0)b'(e),d(e)) = (P0,'0)d'(e),d(e)) and {(f>,v ■0)d'(e),d(e)) = -(« • '0)d'(e),d(e))- 

All the methods that were developed for Dirac spinors can be used here, provided that we assume 
that we are working with E- valued distributions. The structure of the fundamental solutions for 
the wave equation are the same: 

G±(p) = 4{p)Tp{q)S{r,{p)) + V,{p)H^{p). 
where is a bidistribution in £(E) H D'(E) which satisfies the wave equation: 

ip'rGf{p)=M<i): 

Sp{q) being the Dirac mass in p. The application Tp{q) satisfies the equation: 

(Tp(p), 0) = (f> and VjTp{q)V\{q) = 0. (3.3) 



Remark 3.4 The functions r and V are more complex to write and we do not even try to do so, 
since the properties given by the equations I13.3\} are sufficient to conclude. 

A direct consequence of the previous remark is the following proposition: 

Proposition 3.5 The Dirac operator p acting on sections of the fibre bundle E admits two funda- 
mental solutions Gp{q) — p'G^ (q), in D'(E) MD'(E), with, respectively, support in C^{p), for any 
given p, which satisfy, in the sense of distributions: 

p'G^{q)^-5p{l)- 

Finally, we present the compacted spin coefficient formalism introduced by Penrose and Rindler 
m [H]. Let 0-^,^-4 be a given normalized spinor basis and consider the rescaling, for A in C: 

o^^Ao^,i^K-> ^. (3.4) 



Definition 3.6 A spinor <f) is said to be of weight [p, q) if and only if under the transformation 
p.4p . (j) is rescaled as: 

4> I — > yVcf) 

The integer \ {p — q) is the spin-weight of (p and ^{p + q) is its boost-weight. 

We consider the Newman-Penrose tetrad {l,n,m,rn) associated with o^,l^. We define the differ- 
ential operators with regard to these weights: let be a {p, q) spinor. We define: 

p0 \7 ~ pec/) ~ qecf) 

where e — l^\7ioa and a — L^^S/mOA- 

Though the formalism of the Newman-Penrose tetrad will still be used, the usual notations 
o^, t"^ for the basis spin basis are replaced by e^,e^. All the calculations will be performed using 
these notations. We must recall what is the link between these two notations: the spinor basis 
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(o'^j L"^) is rewritten {eQ,ef), so that the dual basis is (e^4,£^) with ~ ~la et e\ — o^. In this 
formahsm, the spinors £^ satisfy: 

J A _ rJ 



Let now consider the field equation for spin 



2 ■ 



V'^^>AB...F = 



for a symmetric field 4'ab...f = 4>{Ab...f) with n indices; for j in {0, 1, . . . , n}, we define: 

— Eq . . ^.EQ^ef . . ei^ <f)AB...F 

n-j times j times 
n-j times j times 

which are the only relevant components to calculate the field (t>AB...F wich can be written, because 
of its symmetry: 



n-j times j times 

= 0j t(A ■ • • (.C Od . . . OF) 

^ n-j times j times 

The quantity (f)j is a (n — 2r, 0) scalar field. It is known to satisfy the following lemma (see [21], 
4.12.42): 

Lemma 3.7 Let j be an integer in {2, . . . , n — 1}. 

Then 4>j+i, , (pj-i and (j)j-2 satisfy the following relation: 

V(f>j - 9'0j-i = {j - l)crVj-2 - j-rVj-i + in- j - l)p4>j - (n - j)K(f)j+i. 

Remark 3.8 This is the more accurate way to write down the constraints equations on the cone, 
since the restriction to the tangential derivatives is obvious. 

We conclude this section by giving the following relation between weighted scalars and differential 
forms (see [25, 4.14.70): 

Proposition 3.9 Let Y. be a two dimensional spacelike closed surface with volume form /is and a 

a (1,-1) weighted spinor. 

Then the integral of S'a over S vanishes: 



j g'a/is = 



3.2 Integral formula for spin | 

Let us consider the future characteristic Cauchy problem for the Dirac operator on E: 

{ pu = Oonr(po) 

\u ^ 0one+ipo) ' ^ ' 

where is a smooth compactly supported function on the cone C~^{pq). It must be noted that the 
problem (|3.5p . contrary to the problem stated in p.ip . does not contain symmetry assumption. 
This assumption will be made afterwards to obtain the integral formula for p.ip . 

By doing the same calculation as for proposition (|2.2|) , a direct consequence of proposition ()3.5|) 
is the following integral formula: 
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Proposition 3.10 Let u be a solution of iS.5\) in E. Then u can he written: 



D(g) 

The formula must now be simplified using the previous methods and a decomposition of the 
spinor u on the same basis as in subsection l2.2.1l u can be written: 

u = (f)A...FeA ■■■ef + V'^F F^A'ci ■■■£f- 

The solution of 



^^^'uAB...F = Oona+(po) 
UAB...F = Oab...f on e+(po) ' 

obtained by projecting on Sa...f the integral formula given in proposition 13.101 
Proposition 3.11 Let ua...f be a solution of: 

r V^^'uAB...F = on a+bo) 

1 UAB...F = Oab...f on e+(po) 

Then, ua...f can be written: 

UA...F - / (hM\ (vf0oB...f(p) -2p0oB...F(p))(VV^o)e^'(9)ei(9)...4('7)/^. 

Ja{q) XT/ 

^{p)VaA' ((^) sfiq)s^s{q) . . . 4(g)) 



(3.6) 



(9) 



(PQB...F[P)y AA' 



A*<T(g)(p) 



+ / V«ro • (V^To • u, V,iq)) ^^^f^ + / (V^To • P'Vq)^iro {p) 

J<y(q) JVj, 

Remark 3.12 Since our interest is in the singular part of the integral representation of the solution, 
we do not give a more explicit expression of the smooth part of the integral formula. 

Proof The first step is to calculate the contraction V^Fq • u: 

V^To • u = -2iV2roe^4' (0a...f£a . • . 4) + 2^V2ro4£!\, (V-^B.^Fel'sl ■ • ■ 4) 
= 2iV2roi-e^ef (j)A...Fe'X + £A£\'^^B...F^i')^B ■ ■ - 4 
= 2iV2rn{(f>oB...F£o' " -b...f^a)£b ■■■^f- 

Since Tp{q) is obtained by doing a tensor product between an element of the spin basis at a point 
p with the spinor obtained by parallely transporting this spinor along the geodesic from p to q, 
which is an element of the spin basis at q, the symplectic product (Tp(g), VTo ■ u) realizes a switch 
between the variables p and q: 

(Tpiq), yPTo ■ u) = zV2ro(0oB...F(p)eo '(<?) ^ B...F(p)£A('?))£i('Z) ■ ■ • 4i<l)- 
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The interversion of the symbols J and p gives (we only make the calculation on Sa...f)- 

iV2\/%,l [ 2iV2ro(f>oB...F{p)4'iQ)^BiQ)---4{'l)t^r„.rA (3.7) 

\Ja{q) J 

= - [ 0oB...F(p)Vi^, f^£jf'(g)ei(g)...4(<z))A*.(,) 



^(Vf0oB...F(p))(VV^o)£jf'(g)£i(g) . . . 4(g)A*.(,) 



which can be simplified in: 



(3.8) 
(3.9) 

T (^) (0oB...F(p))(VV^o)e^'(9)ei('z)...4(9)M.(,) (3.10) 
+ / 2p^(0oB...F(p))(V^^,ro)£^'(g)ei(g)...4(g)A.,(,) 



V 



iV2VV ( / 2zV2ro(^oB...F(p)£^'(9)£i(9) • . . 4(g)/xro.r, 

\J<y(q) 



(3.11) 



(3.12) 



/ (Vf(/)oB...F(p)-2p0oB...F(p))(VVo)£o'('?)ei('z)-..4(g)M.(,) (3.13) 



I^^^c^ob.MpWaa' (^(^^)ef (9)ei(9)...4('?)) M?) 



(3.14) 



The next part of the integral formula is exactly the same as in the case of the Weyl-Dirac spinors 
and is obtained in a similar way.(^ 

Finally, to obtain a solution of the full problem with symmetry, it is sufficient to symmetrize 
the unprimed indices in the formula; we then give a representation, when the problem (13.11) makes 
sense (i.e with adequat restrictions on the curvature for spin greater than §): 

Theorem 3.13 Let ua...f be a solution of the symmetrized characterictic Cauchy problem 

uq = 9ab...f on e+(po) 

where uab...f satisfies the symmetry conditions: uab...f — U(^ab...f)- 

Then the singular part of the integral representation of ua...f, that is to say the part supported on 
the intersection of the cone is given by the formula: 



(Vf<^0B...f(p) - 2p0oB...f(p))(VV^o)£^'(<z)£i('z) . . .4(g)M.(,) 
3.3 Integral formula for spin ^ in the flat case 

This subsection is devoted to the recovery of the Penrose formula; with the same notations as before, 
the following proposition holds: 

Proposition 3.14 (Integral formula for the fiat case for spin ^) Let <j>A...F be a solution of 
p.lSp on the Minkowski space time. 
Then 4> can be written: 

cl>A...F = (-1)" / (Vz(/.o - (n + l)p0o)iA . . . iF^ 
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Remark 3.15 The formula which is given here agrees with the one obtained by Penrose in 
(formula 4-9). The (—1)" comes from the fact that Penrose chooses the convention: 



because of the different choice of normalization (formula (4-7), op. citj; 



whereas our convention is: 



LAO"^ = 1 



OAi"^ = 1. 



Proof : We summarize the geometric elements required to perform the calculation: 
Remark 3.16 We recall the main properties of the spinor basis which was constructed in section 

a. the spinor s o'^ and are constant along a generator of the cone tJ^(po); so that the spin 
coefficients corresponding to the derivatives of o^^i^ along the vector 1°" — o^o"^, K,e,T' are 
zero; 

b. furthermore, for q in d~^{po), the basis {o^, l^) is parallely transported along the integral curves 
of L^T^ and so, in the flat case, is constant along the null generators of the cone 8~ (q); 

c. the derivatives along m of o"^ and are calculated explicitly (see \21f . 4-l<^-l^8): 

g'o^ = -pi^ and^'i^ = -a'o^- 

d. the derivatives of and r can be explicitly calculated by differentiating the relation: 

Pop" = t-q/" + ri-^I^' 
for any p in 3^{po)- Their derivatives are: 

BB'i^^ = --lbob'O^ and \/ aa'T = oaoa' (3.16) 
r 

and, consequently, the only non-vanishing derivative of is 

r 

and the spin coefficients 

t' — —L^ViLA, f' = —L.'^VrnLA, P' = " Ck = —i^'VjnLA and (3 = —a' = —i^^Vml^A 

vanish. 

e. Using equations (|3.16|) and since is a {—l,0)-spinor and r is a (1, 1) scalar, the following 
derivatives vanish: 

g't^ = and gv = 0. 

For the sake of clarity, the calculation is first performed for the Maxwell equations and then for 
the arbitrary spin. The first step is to write the Maxwell equations 

V^^Vab = 
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as 

- Vm(t>a = (tt - 2a)(/)o + 2^01 - k4>2 

Vi02 - Vm-'/'i = -A(?!)o + 2TT(f>i + {p- 2e)(f)2 ,„ 

^ni4>l - V„(/>0 = - 27)00 + 2x01 - (T(f)2 

Vm02 - V„(/ii = -i^(/>o + 2^01 + (r - 2/3)4)2 
with the convention ^oo = 4>o, 4>io = 0i and — 02- We then consider the singular part: 



r 



k (n)\ I 
— ' ''Vf 0ob(p) - 2/90obb)) (MM')eo «^BM<T(g) 



(9) 



j (Vf 0ob(p) - 2p0ob(p))MesAi<T(9) 



with K — Using the first MaxweU equation (|3.17p . and since, for the choice of basis which was 
previously done, the spin coefficients k — C^Vjo^ and tt — —v^'^ila vanish, we obtain: 

^ ^ I { ) (~^) ~ 2P'/'0l(p))'-AOBMo-(g) = J ^ ^ 2^^^^*^™ ~ 2a0oo) MOB Ai<T(?) 

Vwr — r - 2a— 1- Vjnr— — 5 fia{q) 

(•^•j \ 2'Kr J 27rr 27rr^ 2nr 27rr ^ ' 

Since 

V^r = t^o^', V^s,r = t^o^'oBOB' and V^ob = -p^B, 

and 

yp ( 000 M Ob ^ _ 0OOMOB 



by Stoke's theorem (cf. (4.14.70) in [21 ; it is possible to reinterpret this expression using the 
compacted spin coefficient formalism), we obtain: 



B 



<Poo 
(9) 27rr 



We finally have the expected integral formula for the Maxwell equation 



Vf^obb) - 2/9<?!)ob(p))(V^^,ro)eo ('?)£b(9)m^(«) 

{Vicj)oo-2p)iALB^-B 
Zirr 

(Vf0oO-3p).AiB^. 



<t(9) 



-(9) 

The first step of the general proof is to notice that, as in the Maxwell case, the only remaining 
term in the flat case is the equation ()3.13|) since the term p.l4|) vanishes. So the simplification of 
the equation l|3.13p can be done using the same methods. 
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A direct consequence of remark 13.161 is that the relation given in lemma 13.71 is considerably 
simplified: 

Vj e {1, n - - 2p4>, = gV,_i + (n - J - l)p4>,. (3.18) 

In the Minkowski case, the only non-vanishing term in the integral formula is the following: 



which can be simplified in a flat space as: 



Vf0ob...f(p) - 2p0ob...f(p))(V^^^,ro)£^ (9)es(9) . • . ('7)a*<t(9) 



- / TJ^(Vf(/)0B...F(p)-2p(^0F...F(p))i(A(9)£i(9)...£F)(9)M^(9) 

Consider the generic term in this sum: let j be an integer in{l,...,n — 1}: 

/( ) 2^'-^''^^' ~ 2^)'^(^(«)^s('?) • ■ • ^%){(l)^<y[q) 

Since is obtained by contracting j times (f>A...F with l"^ and n — j times with o"^, it means that 
there is exactly j times oa and n — j — 1 times —la in the list e^{q) . . .e^{q); since the sum is 
symmetric, it can be written: 

n~j terms j terms 

Using equation p.lSp . it becomes: 



/ (Vl^j - 2p(l)j)L(^A ■ ■ ■ l-COD ■ ■ ■ OF)ficr{q) 

I 77— - j ~ '^)P<l>jHA ■ ■ ■ ''COD ■ ■ ■ OF)i^a{q) + / -^S'((/)r-l)tA . . . I^COd ■ ■ ■ OF)fJ' 



la{q) ^-"l Ja(q) 

Using remarks 13.161 the last integral is written as a difference 



— O [<pj-l)t(A- ■ ■I'COD ■ ■ ■OF)tJ-a{q) = J O \^ — J f^<T{q) (3.19) 




a{q) 27rr 

+j I p<l)j^lL{A- ■ ■ l-DOE ■ ■ -t-F) ^l^- (3.20) 
Jaiq) 

It has already been noted that: 

a. r is (1, 1) scalar; 

b. (pj-i is a (n — 2j + 2, 0) scalar; 

c. if^A ■ ■ ■ '-cod ■ ■ ■ Of) is a {2j — n, 0) spinor. 

As a consequence, the term integrated in the left-hand side of equation (|3.19p and under the deriva- 
tion S' is (1,-1) spinor. In order to apply lemma [3?9l this spinor is contracted with n constant 
arbitrary spinors; this gives: 




27rr 

n—j terms j terms 
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Finally, we obtain: 

/ (Vi(/>j - 2p^r)i(A. ..(-COD • --Of) = (3.21) 
Jaiq) . , > 27rr 

n—j terms j terms 

[n-J-l) p(t>jl'(A---l'COD---OF)-r^+J p^j-iL(^A---''DOE.-.iF) 
Ja(q) ^ ^ ' 27rr J„r ) s ^ ,s ^ . 



j(q) ^ , ^ y ^Trr ^. ^ ^ 2-Kr 

n-j j n~j+l j-l 

Theses terms are added to obtain the complete expression of the integral formula: 

(V;</.0b. . . f - '2py^A{q)el {q) ■ ■ ■ 4) (?) ^ 

~J \ 3 J Ja(q) ^ V V ' 

n—j terms j terms 

{Wicl>o-2pcl>o)iA...iF^+ 
a(q) 27rr 

X:(-l)"-^-^(n - i - 1) f . P'i'^ HA---^coo... OF) ""^'^ 



n-1 



~( \ J J Ja(q) , V- 



Lf ^'"^'^^ 

,•=1 X . / .<T(g) ^ 27rr ■ 

n-j-l j-1 



The sum is split in two and reindexed: 

n-l 



Y.{-ir-=-\n-j-i)r . M / p4>,.^A...icoo...oF)^''^'^ 



J J Jcr(q) 

n-j j 

n-l 



^■=1 \ 3 J Ja{q) ^ V V ' ^'^'^ 



(3.22) 



The remaining terms are then: 



-^(9) V ' '^^^ 

^0 ""^ 



(-1)"-! / (V,0o-2p<^o)M....^^-(-l)"-lf"^ / p^Po^A-.-OF^^^'^ 



a{q) 27rr V 1 / 27rr 

and the integral formula is, because of the antisymmetry of the symplectic product: 



f)A...F = (-1)" / (V(0o - (n + l)p0o)M . . . (3.23) 



is proved.® 
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Concluding remarks 

a. Klainerman-Rodnianski state in ^14j that a metric (or a square-integrable Riemann curva- 
ture) is sufficient to write the singular part of the Kirchoff-Sobolev parametrix for the Einstein 
equations. It is expected that such a regularity will not prevent the use of this method for 
the arbitrary spin Dirac equation. 

b. The construction of the representation formula is flexible enough to be used with other fiber 
bundles. Provided that the correct geometric hypotheses are stated for the manifold, such a 
representation can thus be obtained for the Rarita-Schwinger (or gravitino) equations. 

c. Chrusciel-Shatah obtained in ,2 L^-estiniates for the Yang-Mills equations. We hope that 
such estimates can be obtained using the Friedlander construction of an integral formula. 
Nonetheless, it must be noted that they intensively used the gauge freedom which exists for 
the Yang-Mills equation: they used both the Cronstrom gauge (to obtain pointwise estimates) 
and the temporal gauge (to obtain estimates on spacelike slices). Similar estimates for the 
Dirac equations could help to explain, for instance, the loss of regularity observed in the 
characteristic Cauchy problem in [Tlj (section 6). 

This work was partially supported by the ANR project JC0546063 "Equations hyperboliques 
dans les espaces-temps de la relativite generale : Diffusion et resonances." 
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